Dissipative Dynamics of Matter Wave Soliton in Nonlinear Optical Lattice 
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Dynamics and stability of solitons in two-dimensional (2D) Bose- Einstein condensates (BEC), with 
low-dimensional (ID) conservative plus dissipative nonlinear optical lattices are investigated. In the 
case of focusing media (with attractive atomic systems) the collapse of the wave packet is arrested 
by the dissipative periodic nonlinearity. The adiabatic variation of the background scattering length 
leads to metastable matter-wave solitons. When the atom feeding mechanism is used, a dissipative 
soliton can exist in focusing 2D media with ID periodic nonlinearity. In the defocusing media 
(repulsive BEC case) with harmonic trap in one dimension and one dimensional nonlinear optical 
lattice in other direction, the stable soliton can exist. This prediction of variational approach is 
confirmed by the full numerical simulation of 2D Gross-Pitaevskii equation. 

PACS numbers: 03.75.Lm;03.75.-b;05.30. Jp 



I. INTRODUCTION 

The dynamics of optical and matter wave solitons with 
different type of management of system parameters has 
been under intensive investigations in the last years 
Two types of modulations have been considered: disper- 
sion and nonlinearity management, which can both oc- 
cur in time and space. Temporal strong and rapid mod- 
ulations of the dispersion are more interesting in opti- 
cal fibers due to many advantages of dispersion managed 
solitons for optical communications [1,13,01 and storage of 
information. Temporal modulations of nonlinearity are 
promising in fiber ring lasers and Bose-Einstein conden- 
sates (BEC) @, 0, i- In the latter case, the suppres- 
sion of collapse, implying in the existence of stable mul- 
tidimensional solitons in attractive BEC, and the gen- 
eration of periodic patterns of matter waves, have been 
predicted @, El, [H, [11 [3 [S [B In optics, 
nonlinearity managed solitons have also been observed, 
as described in Refs. [H, [H, [H, [2l| . 

Recently, the attention has been devoted also to the 
periodic spatial management in nonlinear optics and 
Bose-Einstein condensates. In optical media, the non- 
linear Kerr coefficient can be periodically modulated in 
space, leading to the problem of optical beam in a 2D 
medium with nonlinearity management. In BEC, the 
spatial variation of scattering length is possible [22|, [H, 
[13, [H [H [13, [H, for example^y using optically in- 
duced Feshbach resonance [2^, [30] ■ In elongated con- 
densates new types of localized matter waves packets 
can exist. In 2D case, the situation is less clear. The 
study of one-dimensional (ID) nonlinear periodic po- 
tential in two-dimensional (2D) non-linear Schrodinger 
equation (NLSE) shows that broad solitons are unstable. 
As verified in Ref. [31], narrow solitons centered on the 
maximum of the lattice potential can be stable, but the 
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stability region is so narrow that they are physically un- 
stable. Stable gap solitons can exist in the BEC under 
combination of a linear and nonlinear periodic potentials 
[32, 33, 34]. 

However, models considered till now are strongly ide- 
alized. In particular, using the optically induced Fes- 
hbach resonances we can generate mixture of conserva- 
tive and dissipative nonlinear optical lattice. In view of 
that, around the Feshbach resonance, one can observed 
non- vanishing contributions to the imaginary part of the 
scattering length. 

In this work, after an analysis of a conservative system 
with nonlinear optical lattice, we will consider the influ- 
ence of nonlinear dissipation on the dynamics and the 
stability of solitons. In particular, we note that the role 
of such kind of dissipation can be crucial for the existence 
of solitons in multi dimensional nonlinear optical lattices. 
Such hope is supported by the well known fact that the 
homogeneous nonlinear dissipation can arrest collapse in 
the cubic focusing multi-dimensional NLSE [3^. The 
possibility of existence of dissipative solitons is investi- 
gated, considering compression effects and atom feeding. 

The organization of the paper is as follows. The model 
is described in the next section. In Sec. 3, it is investi- 
gated the properties of localized states in case of attrac- 
tive and repulsive 2D condensates in ID nonlinear optical 
lattice, with or without harmonic trap in one of the di- 
mensions. In Sec. 4, it is performed an analysis of the 
evolution of 2D soliton under ID periodic nonlinearity 
and dissipation, using the variational approach and by 
direct numerical simulation of the GP equation. 



II. THE MODEL 

Recently, it has been suggested to generate nonlinear 
optical lattices in BEC by two counter propagating laser 
beams near the optical induced Feshbach resonance (25l . 
[26j . The spatial variation of the optical intensity leads 
to a spatial periodic variation of the atomic scattering 
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length. Such structure can support new types of localized 
nonlinear states. The GP equation for the wave function 
■0 = 'il){xi, X2,t) has the form 



2m 



dxl 



5(2:1, a;2)|V'P'/', (1) 



dt 
where 

g{xi,X2) = Y + (gi + 152) (cos^(fcxi) + cos^(fca;2)) . 

(2) 

50 is related to the s— wave two-body scattering length 
tts, go = (47rft^/m)os, with go > (go < 0) for attractive 
(repulsive) condensates; gi (> 0) is related to the optical 
intensity; and 52 parametrize dissipative effects. 

The optically induced scattering length and the inelas- 
tic collision rate coefficient i^ine; (imaginary part of a^) 
are described by [29| 
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where Uso is the scattering length without light, A is the 
detuning from the photo-associated resonance, and hki 
is the relative momentum of the collision. TsUm is the 
resonant transition rate between the continuum state and 
the molecular state, proportional to the laser intensity 
I{x). Far from the resonance, the imaginary part of the 
scattering length is small, such that 3(a) ^ di{a). It 
was shown in Ref. (s^l that, in the experiment with ^^Rb 
one can obtain optically induced large variations of the 
scattering length. The laser intensity was 460 W/cm^ 
and the variations occurred from lOao to 190ao (with 
Oso = lOOao and oq the Bohr radius). 

By considering the following variable changes and def- 
initions in II]), 



Kx — 2kxi, Ky — 2kx2, t — — 



er ^ h^e 



7i=0,l,2 



2|5o|' 



we obtain the dimensionless equation 



.du d^u ^d'^u , ,2 



(5) 
(6) 

(7) 



where 



7(0;,?;) EE 7o-(-(7i+i72)(l+(5o+cos(Kx)-|-^ocos(Ky)), (8) 
and the wave-function was redefined such that 



u = u{x,y,T) 



(9) 



From ([61) to ([9]), we should note that 70 is fixed to —1/2 
for attractive condensates; and -1-1/2 for repulsive con- 
densates. Different cases can be considered: ID geom- 
etry is realized when /3 = Jo = 0. Anisotropic 2D case 



is realized for /3 = 1, = 0. And the 2D isotropic 
case can be achieved with /3 = = 1. Next, we consider 
more explicitly in our study the anisotropic 2D case, with 
—\ and — As the soliton is completely free in the 
y— direction, we also examine the possibility to have a 
harmonic trap mcL;|a:|/2. Following the transformations 
(III), a dimensionless frequency lo is also defined, such that 
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In order to extend our study of the stability conditions 
in a few realistic cases, it is also verified the effect of 
a compression, which can be achieved by an adiabatic 
time variation of the background value of a scattering 
length [22], by modifying 70 as 



70 70 (t) = 70 exp [2a (r - r^)] 6'(r - r^). 



(11) 



Compression effect can also be achieved by a feeding pro- 
cess, which can be described by an additional term iafU 
in the GP equation [3^. Note: If the modulation of non- 
linearity in time is induced by increasing the transverse 
frequency of the trap, then we should multiply the full 
nonlinear term by exp(2Q!r). With these considerations, 
the Eq. ^ can be written as 

.du d^U d^U , ,2 2 2- ^ir.\ 

^TT = -^-^-7(a^,y,T)|'«| V u+iafU, (12) 



where 



dx^ dy 



y, t) = 7o(t) + (71 + 172) [1 -I- cos(kx)] 



(13) 



In the above, one should take a/ = when a 7^ in (jlip . 
as such parameters have similar role in the formalism. 



III. CONSERVATIVE SYSTEM 

It is useful to describe shortly the solitons and their 
stability in the conservative case (72 = 0). One dimen- 
sional conservative case has been considered by using a 
variational approach (VA) in [2^. Using an exact ap- 
proach, the 2D case with ID nonlinear optical lattice was 
studied in jsij, where the authors have considered the 
case with attractive background nonlinearity (71 > 0). 
Looking for perspective applications to BEG, we will con- 
sider here the 2D problem with 1 D nonlinear optical lat- 
tice. Following Ref. [26], we start our analysis using the 
VA formalism. 

With m(x, y, t) = v{x, y) exp(— i/^r) in (fT2)l . and taking 
a/, a and 72 equal zero, we obtain 

cP'v cP'v 

^ ^ = (70 +7iCOs(Ka;))u^-|-w^7;\, (14) 

where we are redefining 70 to 70 = 70 + 71 . In view of 
our definitions in ([51), this implies that for attractive 
condensed systems we have 70 = 71 + 1/2; and, for 
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repulsive ones, 70 = 71 — 1/2. The sign of 70 gives the 
sign of the background field. But, we should note that 
we can have situations where the same 70 > can refer 
to attractive or repulsive condensates. Example: 
7o = 1 with 71 = 1/2 and 70 = 1/2 (attractive); 
70 = 1 with 71 = 3/2 and 70 = —1/2 (repulsive). 
These two situations differ in (fH)) . because the strength 
of the oscillatory term is different. However, as the 
results are similar, we prefer to analyze separately the 
cases of repulsive condensates with negative background 
field, which occur when < 71 < 1/2 (0 > 70 > -1/2). 

The corresponding averaged Lagrangian L is obtained 
from the density £, as 



L 



— oc 'J —00 



dx / dy 



(15) 
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dx 




dy 



70 + 71 cos(fi:a:) ^4 



Here, it is interesting to observe that a scaling given by 
K is applied to the observables obtained from the above 
equations, as the root-mean-square radius in x and y di- 
rections, chemical potentials, frequencies and energies. 
In order to see that, we can redefine all the observables 
using the variable transformation, x = nx and y = ny, 
such that we have no k dependence in a new set of ob- 
servables (represented with a "bar" ) that are being calcu- 
lated. This scaling essentially implies to consider k = 1 
in all the equations. At the end, the physical observ- 
ables will be given by the relations ([5]) (with k = 1). For 
example, in the case of mean-square radius we will have 

4fc2 



and 



= (17) 



Next, in our VA we consider the Gaussian ansatz 

.(x,y)=Aexp(-^-|^ 

which has the normalization given by = 770102^42 
corresponding averaged Lagrangian is given by 



(18) 
The 
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From the Euler-Lagrange equations for the parameters, 
dL/dN = and dL/dai=i^2 = 0, we obtain: 
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+ w^a^, (20) 
(21) 

- 1 = 0.(22) 



In the case that w = 0, this set of equations, for /i, N ^ 
and 02, can be expressed in terms of ai, as 
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(23) 



(24) 



435) 



For the case that cj 7^ 0, the relation for in terms of 
fli can be obtained from (|22p and 



a-i 



cja2,o 



(26) 



Equations P^ . and (|2ip form the set of equations 
for 7^ 0. 

Next, we consider separately the cases of attractive 
systems, with 70 = l\ + \ > 0, and repulsive ones with 

70 = 71 - 5 < 0- 



1. Attractive condensate (70 = 71 + 1/2) 

This case, which corresponds to 70 = 1/2 and 71 > 0, 
has been investigated recently in 31]. With w = 0, it is 
applied the set of equations and P^ . In the 

general case with w 7^ 0, we should consider Eqs. (|20p . 
([21]), and (IMl). 

Let us first verify the analytic limiting cases of the VA 
expressions, for cj = 0: 

12,0 ~^ for ai << 1 and a\ >> 1; 

— 1/a^, for ai << 1 and a\ » 1; 
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for ai = 0; 



70 + 71 271 + 1/2 

AT ^TT An 

70 71 + 1/2 
Limiting cases of the VA expressions, for w 7^ 0: 



a2 
A^ 



ai, for ai << 1; 
1/y/uj, for fli >> 1 

— l/af, for ai << 1 

— 2/0^ + — > for fli >> 1; 

for fli = 0; 



70 + 71 ' 



AT '1^"2 „ r 

N > 0, tor oi 00. 

7oai 



In Fig. [U we plot the corresponding results for the 
chemical potential /i as a function of A^ (upper frame) and 
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iV as a function of ai (lower frame) . Numerical soutions 
to PDE results were done with algorithm presented in 
Ref. [33. Considering the Vakhitov-Kolokolov (VK) cri- 
terion [38| for the soliton stability, dfi/dN < 0, from the 
results given in the upper frame of Fig. [Tl we note that 
the solitons are unstable. This result is in agreement with 
the prediction of Ref. [sT'l. We note, from the VA results, 
that in the limit of large ai the system has a tendency to 
stabilize, indicating that with just a small trapping po- 
tential we can produce a stable region. This behavior is 
shown by the VA results given in Fig. [21 The variational 
approach, besides an expected small quantitative shift, 
provides a good qualitative picture of the results when 
compared with full numerical predictions. If one is first 
concerned with the stability of the system (instead of the 
quantitative results of the observables) , the VA provides 
a nice and reliable picture. 
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FIG. 1: Attractive case, with 70 = 1 and 71 = 0.5. Re- 
sults for the chemical potential /i, as a function of A'^ (upper 
frame) and A'^ versus ai (lower frame), obtained using vari- 
ational approach (VA) and full numerical calculations. The 
variational parameter for the width, ai, and the root-mean- 
square radius, 



(2;^), are related by ai — ^J2{x'^) [Actually, 
the physical observables depend on k as given by ([5} and (|17j . 

In our VA, when we keep uj fixed (zero or nonzero) 
and increase the value of 71, we observe that the gen- 
eral picture in respect to stability of the system does not 
change. This lead us to conclude that we cannot improve 
the stability of the system by increasing the strength of 
the lattice periodicity for attractive condensates. In the 
following, we are going to analyze the cases with 70 < 0. 
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FIG. 2: Attractive case, with 70 = 1 and 71 = 0.5. VA results 
for the chemical potential /i, as a function of N (upper frames) 
and N versus a\ (lower frame). The results are given for /i 
near zero, considering three values of the frequency; a; = 
(solid line), 0.01 (dashed line) and 0.045 (dotted line). 



2. Repulsive condensate, with 70 < (70 = 71 — 1/2). 

We should remind that by repulsive condensate we 
mean an atomic system where the particles have origi- 
nally positive two-body scattering length, such that in 
(2]) we have < 0; or 70 = —1/2. So, given 71 (pa- 
rameter of the spatial periodic variation of the atomic 
scattering length), 70 = 71 — 1/2. And, if we also con- 
sider a negative background, such that 70 < 0, 71 will be 
restricted to the interval < 71 < 1/2. 

Some other limitations are applied in the parameters, 
considering that the widths and N must be real positive 
quantities. The relation between the widths 02 and oi, 
Eq. ((24| for = 0, implies in a limitation to the values 
of ai: 



> 



1 



-1, 



7^ = — In 

l,rnax o 



71 



71 



(27) 



This limit, ai^rnax, is necessary in order to have 02 and 
N real and positive quantities for any values of w. It will 
also restrict the actual values of the parameter 71 to 1 /4< 
7i <l/2. The cases with 71 > 1/2 are also allowed, in 
principle, without upper limit for ai. However, such cases 
will correspond to positive background field, 70 > 0, that 
have already been considered in the previous subsection. 



5 



In view of the above, let us also verify in this case the 
analytic VA limits. For cu = 0: 



a2,o 



Mo 



And for cj 7^ 



ai, for ai << 1; 
oo, for ai = ai^max- 

-1/al 
l/{2al), 

r 8^/(471 - 1), for ai = 0; 
1^ oo, for ai = 



for ai << 1; 
for ai = ai m 



a2 



N 



0.1, 

l/(2a?) 



for 
for 



fli = 0; 

fll = CL\ rnax- 



fli << 1; 

Q-l — ai^rnax- 

for fli << 1; 
for ai = ai^maa;. 

87r/(47i - 1), for 
327r/[(l-27i)VJJK2a?^,„„J, for 

In Fig. [21 we plot N versus ai and the chemical poten- 
tial fi versus N, for 70 = —0.1 and 71 — 0.4, considering 
VA and four values of uj (0, 0.07, 0.1, 0.3). In the case 
of w — 0, we also include results obtained from exact 
PDE calculations. Following the VK criterion for sta- 
bility, d^/dN < 0, we notice that stable regions start 
to appear with lo k, 0.1. With w > 0.3k^ the unstable 
regions almost disappear. However, as one can observe 
in the lower frame, the width ai is quite limited due to 
the condition (j27p . The observables and ai depend on 
the wave parameter k of to the spatial periodic variation 
of the atomic scattering length through the scaling rela- 
tions ([5]) and ^ with k — 1. However, contrary to some 
discussions and conclusions of Ref . [3l| , specific values of 
the parameter k cannot affect the conclusions on stabil- 
ity. In such cases of conservative systems, the stability 
results from combined effects given by the parameters 
7o, 71 and w. Our main conclusion is that, without the 
trapping potential (included in the y— direction), taking 
bj — the optical lattice cannot stabilize the solutions, 
neither for repulsive nor for attractive condensates. 

In order to further check the role of the optical lattice, 
for the repulsive case we also investigate the case with 
constant lo and different values of 71. From the results 
shown in Fig. [3l for 71 = 0.4, we found appropriate to 
consider u) = 0.07, which has a marginal stability near 
/J, R:! 0.05. The results are shown in Fig. [H where we first 
observe that a larger 71 can help to allow the width oi 
to increase, within the limiting condition (|27p . However, 
the marginal stability remains for corresponding different 
values of the chemical potential. In order to keep the 
plots of Figdlfor different values of 71 in the same frames, 
we have normalize the number N such that it is equal to 
one when oi is zero. 

The plots of the evolution of profiles are presented for 
LO = Q and 0.3 in Fig. [5l confirming the VK prediction. 
The results using the VA have good agreement with PDE 
prediction, as shown for uj =0 and 0.3. The profile at 
t = 50 is practically undistinguishable from the initial 
soliton form. 
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FIG. 3: Repulsive case, with 70 = —0.1 and 71 = 0.4, for 
fj, versus A'^ (upper frame), and TV versus ai (lower frame). 
In both frames, we show the results using the variational ap- 
proach, for Lj = 0, 0.07, 0.1 and 0.3. In the upper frame, the 
exact PDE results are also shown in two cases: a; = (for 
which the system is unstable) and uj = 0.3 (for which the sys- 
tem is stable) . In this last case, near the region where the VA 
presents a small unstable branch (22 < N < 25), our exact 
numerical results are shown only for A'^ > 24. As observed, 
the VA is giving a general picture of the exact solutions. 

IV. EVOLUTION OF 2D SOLITON UNDER ID 
PERIODIC NONLINEARITY AND DISSIPATION 

In this section, we will consider the case we have 72 7^ 
in (|12p . To study the dynamics of a 2D soliton with 
ID periodic nonlinearity and dissipation, we also apply a 
variational approach and full numerical calculations. In 
the Gaussian ansatz (|18|) we should also include param- 
eters related to dissipative effects and initial conditions. 
For a bright soliton, the ansatz can be taken in the fol- 
lowing form: 



Aexp 



2ai 



2al 



X exp [\bi{x - xof + ib2y^ + i0(T)] , (28) 

where bi {i =1,2) are related to dissipative effects, with 
xq and the phase </> related to initial condition. 
The Lagrangian density for Eq. p2)) is 



C{x, t) = ']^{uru* - <«) - \u.A^ - \Uy\' + ^1^/1^(29) 



where 7(2;) is given by (|13p . Next, from the ansatz 
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we obtain the corresponding averaged Lagrangian: 



L = 



dx / dyC 



" a2 

= A aiflT 

2 



1 



afibir + ibf) + — + ai,{b2r + 46^) 



+ ^ + 2(/.,-^(7o + 7iCOs(^xo)e-«'''?/«) (30) 

The equations for the sohton parameters rji = [A, a, b, (j)] 
in the VA are derived from(see for example [sot) 



dL d dL 



dr]i dr dr]i^r 



dx 



dy[R-^ ^ c.c 



(31) 



FIG. 4: VA results for the repulsive case, with 70 — —0.1 and 
uj fixed to 0.07, considering 71 — 0.35, 0.4, 0.45 and 0.5. In 
the upper frame we have ^ versus N/N{ai — 0); and, in the 
lower frame, N/N{ai =0) versus ai. 
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FIG. 5: The soliton profile in the stable region predicted by 
the VK criterion for lj — 0.3, 70 = —0.1, 71 = 0.4, fi = 0.1, 
A « 33 at the time t — 50. Solid line is x-direction, dashed 
line is y-direction. 



where the perturbation term R is 



R 



-172(1 + cos(Kx))|upu + iafU. 



(32) 



Here we are taking into account a linear amplification 
term (a/) describing the atoms feeding. Finally, from 
the above, we obtain the following system of five coupled 
ordinary differential equations (ODE) to be solved for 
the parameters of our variational approach (VA): 
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iV = 7rA'^aia2, with «, j — 1,2 (i 7^ j), we have: 
2a/ AT, 

l + e-«'a?/8cos(Ka;o) + ^ 

70 + 71 l + e-'''"i/*cos(Ka;o) + ^ 
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(35) 
(36) 

(37) 



(38) 
(39) 
(40) 



In the next, we present some of our results, when 
considering periodic nonhnearity with dissipative effects. 
Considering the scahng of observables with k, discussed 
for the conservative systems in section III, which can also 
be verified in the present case, we have the corresponding 
transformation bi bi/n^. 

In Fig. [51 we have results for the of full numerical simu- 
lations (PDE) for the evolution of the matter wave packet 
under combination of the conservative and dissipative 
nonlinear optical lattice in the case of the attractive con- 
densate 7o = 1/2. As we can see the collapse is arrested 
by the dissipative nonlinear optical lattice. The results 
are compared with the prediction of the VA approach 
(ODE). We observe a good agreement of the VA with 
full-numerical calculations. 

We also have investigated the role of a deviation 5 of 
the given norm from the critical norm, in the initial wave 
packet: A A{1 + S). The results of the full numeri- 
cal simulations are presented in Fig. [T] Increasing the 
deviation S from the critical norm, multiple peaks are 
observed, corresponding to revivals of the wave packet 
during the collapse. The number of peaks growths as 
S varies from 0.02 to 0.5. The focusing-defocusing cy- 
cles connect the action of the periodically varying in the 
space with the inelastic three-body interactions. In linear 
conservative optical lattice, with inelastic three-body in- 
teractions, the focusing-defocusing oscillations have been 
studied in Ref. U^. 

The spreading out of the pulse after the collapse is ar- 
rested, observed in Fig. [51 can be compensated, as we 
note previously, by an adiabatic variation of the back- 
ground scattering length, described by a time variation 
of 7o. When we consider a time dependent 70, as given 
in Eq. pT|) . the feeding term parameter af should be 
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FIG. 6: Results for the amplitude as a function of r, us- 
ing variational approach (VA) and full numerical calculations. 
The VA and the full numerical calculations have the value of 
fixed to the same value for r = 0, implying in a small shift 
of A(0), as shown by the results. 

zero, because one can show (with a redefinition of the 
wave- function) that it has a similar effect. In Fig. [51 we 
show our full numerical results confirming the stabiliza- 
tion of the condensate after the collapse was arrested. 
The mechanism of this stabilization was given by an ap- 
propriate tunning of the parameters a and Tc of Eq. (jlip . 

V. CONCLUSION 

Dynamics and stability of matter-wave solitons in the 
mixture of conservative and dissipative nonlinear optical 
lattices are investigated, considering 2D BEC, with ID 
conservative plus dissipative nonlinear optical lattices. 
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FIG. 7: Results for the amplitude as a function of r, showing 
the collapsing behavior for k = Att. The dissipation is fixed 
to 72 — 0.0025 and 5 is increasing from 0.02 to 0.5. 
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FIG. 8: Soliton stabilization via compression effect, with Tc = 
5 and a = 0.018 in Eq. pip . The other parameter are 72 = 
0.006, S = 0.01, and k = Att. 



In the first part of this work, it was analyzed conserva- 
tive systems, with nonhnear optical lattices, for attrac- 
tive and repulsive condensates. It was clarified the role of 
the scales when calculating the observables as the chemi- 
cal potential and the widths. Our conclusion is that, in a 
2D system, a nonlinear periodic lattice in one direction by 
itself cannot give stable solutions, satisfying the VK cri- 
terium 38]. Such periodic lattice in the x— direction can- 
not compensate the collapsing effect which results from 
the other dimension. We verify that stable solutions can 
be obtained by controlling the soliton with a harmonic 
trap in the y— direction. For repulsive condensates, the 
2D stable soliton can exist in the geometry with ID non- 
linear optical lattices in one direction and harmonic trap 
in the other direction. 

In the second part of the work we analyze the dynam- 
ics of the above 2D system, with periodic nonlinearity in 
X— direction and without trap in the direction, when 
we add non-conservative nonlinear optical lattice terms. 
We show that the collapse of the condensate can be ar- 
rested by a dissipative periodic nonlinearity. To study the 
evolution of 2D wavepacket we apply the time-dependent 
variational approach. To compensate the wavepacket 
broadening, the adiabatic time variation of scattering 
length is used. It is shown that the metastable dissipa- 
tive soliton can exist in 2D condensate with ID periodic 
nonlinearity. Analytical predictions are confirmed by the 
numerical simulations of full 2D GP equation. 
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